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An unbinned maximum likelihood analysis of CDMS low-energy data reveals a strong preference 
(5.7(T C.L.) for a model containing an exponential excess of events in the nuclear recoil band, when 
compared to the null hypothesis. We comment on the possible origin of such an excess, estabhshing 
a comparison with anomalies in other dark matter experiments. A recent annual modulation search 
in CDMS data is shown to be insufficiently sensitive to test a dark matter origin for this excess. 
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The CDMS collaboration has recently made public a 
negative search for an annual modulation in lovif-energy 
signals from their cryogenic germanium detectors 
This effect is expected from Weakly Interacting Massive 
Particle (WIMP) interactions with dark matter detector 
targets [3] ■ Observation of this WIMP signature has been 
claimed by the DAMA collaboration vifith high statisti- 
cal significance ^3:], using lovif-background Nal(Tl) scin- 
tillators. The CoGeNT collaboration recently released 
fifteen months of data from underground germanium de- 
tector operation Q . These display a compatible modula- 
tion [J-|6|, albeit with the smaller statistical significance 
that would be expected from a short exposure. 

Fig. 6 in [ij shows, for the first time, detailed in- 
formation from all eight CDMS germanium detectors 
employed in the modulation search and a previous low- 
energy analysis 0- Specifically, it contains the distribu- 
tion of single-interaction events in the ionization energy 
{Ei) vs. recoil energy (Er) plane that can be used to 
identify their origin in nuclear recoils (NR) like those ex- 
pected from WIMP and neutron interactions, or electron 
recoils (ER) like those induced by gamma backgrounds. 

A formal assessment of the possibility that a significant 
WIMP-like low energy NR population might be present 
in CDMS data is of particular interest. The CDMS col- 
laboration has not made public a dedicated study of this 
type, but has put forward arguments Tj indicating that 
the majority of their low-energy events would originate 
in unrelated backgrounds. These arguments have been 
criticized in [8]. Additional objections will be presented 
here, arising from the analysis described below. 

Fig. 1 displays the sum of these data in Ei vs. Er rep- 
resentation for dark matter search runs collected over a 
two year span. The response across detectors to gamma 
and neutron calibration sources varies from device to de- 
vice, as denoted by the individual centroids of the ER 
and NR bands (thin dotted lines). While these varia- 
tions have the effect of somewhat blurring the distinction 
between background sources, they are sufficiently small 
not to preclude an attempt to analyze summed data, in 
an effort to maximize the available statistics. Besides 
the mentioned ER and NR components, two additional 
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FIG. 1: Scatter plot of single-interaction events in all eight 
CDMS detectors, digitized from individual plots in using 
gSdata software [^. Thick solid lines represent the approx- 
imate centroid of the electron recoil (ER) and nuclear recoil 
(NR) response bands, averaged over all detectors, with thin 
dotted lines corresponding to each individual device, aver- 
aged over the small temporal variations observed during a 
two-year span (Fig. 6 in [l|]). A third thick line indicates the 
nominal centroid for zero-charge (ZC) events, defined by the 
condition Ei — (see text for a discussion on deviations). 
A disproportionate fraction of events in the region Ei < 1.5 
keVee, Er > 5 keVnr, come from T3Z5, a detector with a 
markedly degraded energy resolution. The blue boundary is 
the (detector-averaged) region selected by CDMS for an an- 
nual modulation search [l[. 



sources of background are described in : a diffuse sur- 
face event component (SE, Fig. 1) due to gamma and 
beta interactions on detector surfaces, leading to partial 
charge collection, and a population of zero-charge (ZC, 
Fig. 1) events taking place on the edges of the detectors, 
resulting in complete charge loss due to local electric field 
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distortions. 

We employ the RooFiT library [13], part of the ROOT 
software distribution, to perform an unbinned maximum 
likelihood analysis on CDMS data. The purpose is to 
compare the null hypothesis (a model containing ER, 
SE and ZC components only) with the alternative hy- 
pothesis, which includes an additional NR component. 
The ER, NR and ZC components are defined similarly 
by probability density functions (PDFs) of the form 
A{Er) X G{E,,Er), where A{Er) = Ai ■ e^^^ -^-, Ai, As 
are positive-defined free parameters, and G{Ei,Er) — 
^^(Ei-c{E^)f/2a\Er) ^ rpj^^ addition of a constant term 
to A{Er) is seen to be unnecessary. The function C{Er) 
returns the ionization energy of the centroid for the ER, 
ZC and NR bands (Fig. 1). This function is linear for 
ER and ZC (C(£'r)^0 for ZC, see discussion below), and 
quadratic for NR. The function a\Er) = Sl + S2-C{Er), 
with Si, 5*2 positive-defined free parameters, accounts for 
the energy resolution of the detectors, dominated by the 
significantly noisier ionization channel (a factor ~3 larger 
dispersion than along the phonon channel 11, 13] )■ The 
comparatively small smearing of the data from the resolu- 
tion along the phonon (recoil energy) channel is absorbed 
by the model through A{Er). The functional form of 



a (Er) arises from a conventional description 13|, where 



encapsulates all sources of electronic noise and S2 ac- 
counts for the statistics of information carriers (electron- 
hole pairs here). While ER, ZC and NR components 
of the model are allowed independent values of Ai and 
A2 (i.e., are allowed to grow exponentially towards di- 
minishing Er independently of each other), the energy 
resolution parameters 5*1 , 5*2 are common to all compo- 
nents: in other words, the energy resolution is blind to 
the nature of an event generating a pulse of a given Ei. 
Finally, based on the information in [l|, 0, [13] , we expect 
the combined triggering and software cut efficiency to be 
sufficiently close to 100% over the range of Fig. 1 to be 
conservatively neglected. 

Following the description provided in SE events 
are combined with those in ER through a " Crystal Ball" 
PDF [1^, conveniently pre-defined within RooFiT. This 
approach (a Gaussian core as above to account for ER, 
with a power-law low-Ei tail for SE) describes the no- 
ticeable " bleeding" of ER events into the SE population 
due to partial charge collection at surfaces. The dis- 
tribution of ER events in energy is essentially flat for 
the studied range, with a peak of negligible intensity at 
Ei = 6.5 keVee from cosmogenic activation of long-lived 
^^Fe [1, [13] as the only discernible feature. 

The description above is used for the analysis of indi- 
vidual detectors, but an additional allowance must be 
made for an adequate treatment of summed detector 
data: while C{Er) in principle contains no free param- 
eters, the coefficients defining these linear or quadratic 
functions are allowed a modest Gaussian smear defined 
by the variance in their values across detectors. This 
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FIG. 2: Grayscale intensity plot equivalent to Fig. 1, with 
detector T3Z5 removed. Dotted lines represent the best-fit 
(ilcr) boundaries for the ER, NR and ZC bands. The best- 
fit to the electronic noise of the detectors compares well with 
a detector-averaged value in [T^j (303 ± 30 eV vs. 293 eV, re- 
spectively). We employ the same model description as CDMS 
to account for ZC events (see text). 



accounts for detector differences in energy scale calibra- 
tion, energy resolution, exposure, ER background rate, 
and the necessary averaging of the small temporal vari- 
ations observed in their response [J. We apply the re- 
quired penalty to the goodness-of-fit when enabling this 
constrained fit. 

The dispersion along Ei for ZC events in any given de- 
tector is defined by the electronic noise of the ionization 
channel, and as such should be symmetric around Ei = 0. 
The effective ZC centroid can nevertheless acquire a non- 
zero value due to, for instance, a less-than-perfect energy 
calibration along Ei . We conservatively allow for a posi- 
tive ZC centroid (Fig. 2) even when analyzing individual 
detectors, obtaining excellent agreement with CDMS's 
identical modeling of ZC events. Specifically, matching 
within 70 eV both positive (Ei = +0.7 keVee) and neg- 
ative {Ei = —0.26 keVee) 2a ZC boundaries for detec- 
tor T1Z5, the only one for which this information has 
been made available [l^. In spite of this possible small 
ZC shift, the symmetry of positive and negative fluctu- 
ations in detector noise about its baseline is to be ex- 
pected, and well-known to workers in radiation detection 
17| . In the case of CDMS detectors it corresponds to 
a symmetric "zero strobe" peak [3] triggered by a ZC 
event registered along the phonon (recoil energy) channel 
only. Keeping in mind the physical origin of ZC events 
(low-energy electrons traveling through edge regions) we 
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expect, if anything, less dispersion in Ei about zero the 
lower the energy (shorter the electron range) of the inter- 
action. High-energy ZC events of up to 100 keVnr remain 
tightly clustered about Ei'^O 'tJ . The absence of signif- 
icant ZC band distortion or broadening with decreasing 
Er is evidenced in CDMS ^^^Ba calibrations [lit . 

It is worth mentioning at this point that the ioniza- 
tion yield vs. recoil energy plot for detector T1Z5 offered 
in most recent CDMS public presentations displays 
markedly different boundaries for the ZC band (±2cr) and 
NR band (+1.25(7,-0. 5(t). This, together with the chosen 
Ei/Er vs. Er representation, communicates the percep- 
tion that most low-energy events below ~ 5 keVnr belong 
to the ZC population. A more objective representation 
like that of Fig. 2 results in the immediate impression, 
confirmed by our formal analysis, that a large fraction of 



these events belongs to NR. 

Table I contains the results of comparing null and al- 
ternative hypotheses via the standard log-likelihood ratio 
test statistic, for summed data and individual detectors. 
A is the ratio between the likelihood of the null and alter- 
native hypotheses, and Ad.o.f. represents the difference 
in degrees of freedom. The small p-values obtained from 
this statistic [l^ indicate that the model containing a 
NR component is strongly favored, a statement that can 
be maintained with a > 5a C.L. for summed data. The 
best-fit to the integrated NR event rate and NR expo- 
nential decay constant seems compatible across all 
cases studied, once the uncertainties and detector differ- 
ences (exposure, resolution) are taken into account. Fig. 
3 displays the projections of the data and best-fit model 
components along the Ei and Er axes. 



Detector 


Exposure (kg-day) 


-2 ln(A)/Ad.o.f. 


p- value 


Si(eVee) 


NR rate > 3 keVnr (c/kg-day) 


A^" (keVnr-i) 


all 


240.4 


37.8 / 5 


4.1 • lO"' 


303 


0.93 ±0.17 


0.650 ±0.081 


all-T3Z5 


214.6 


45.9 / 5 


9.5 • 10"^ 


266 


0.95 ±0.14 


0.701 ± 0.079 


T1Z2 


43.4 


13.5 / 2 


1.1 • 10"^ 


220 


0.77 ± 0.20 


0.569 ±0.139 


T1Z5 


35.0 


12.5 / 2 


1.9 • 10"^ 


201 


0.50 ±0.17 


0.714 ±0.299 


T2Z3 


28.0 


5.7 / 2 


5.7-10-2 


246 


1.31 ±0.68 


0.659 ± 0.222 


T2Z5 


34.7 


2.7 / 2 


2.6 ■ 10-1 


442 


0.76 ±0.41 


0.745 ± 0.355 


T3Z2 


7.8 


3.9 / 2 


1.4 • 10-1 


333 


3.38 ± 1.24 


0.705 ± 0.215 


T3Z4 


29.6 


12.8 / 2 


1.7 • 10"^ 


142 


0.39 ±0.15 


0.636 ± 0.202 


T3Z5 


25.8 


0.22 / 2 


8.9 • 10-1 


406 


0.15 ±0.38 


2.01 ± 2.44 


T3Z6 


36.1 


3.27 / 2 


1.9 • 10-1 


228 


0.61 ±0.31 


0.707 ±0.361 



TABLE I: Comparison of null and alternative models via the log-likelihood ratio (see text). As expected, the p- value for 
individual detectors correlates to their electronic noise (Si), responsible for loss of energy resolution, and with it the ability to 
distinguish between background components. 



The physical origin of an exponential excess along the 
NR band, concentrated below a few keVnr only, is a 
daunting question. A few remarks can nevertheless be 
made with some degree of certainty. CDMS has re- 
cently stated Q that Monte Carlo simulations produce 
a negligible < 1 event background for the dataset exam- 
ined here. We agree with this statement: based on Co- 
GeNT simulations, scattering of underground neutrons 
from cavern and (/i,n) origins is expected to produce a 
much less steep rise in NR rate towards low Er, and 
a low-energy NR rate negligible when compared to the 
magnitude of this excess [20] . In earlier searches [2l'| , the 
CDMS collaboration found no evidence for NR events 
above 15 keVnr (an aggressive SE cut is only possible 
above ^ 10 keVnr [H, 0]), confirming that a neutron ori- 
gin is highly implausible. Separately, events registering 
in multiple detectors ("multiples". Fig. 6 in [l|) are seen 
to spill into the low-energy NR bands. A fraction of 
these originate from surface events involving simultane- 
ous beta and gamma emission (e.g., from ^i"Pb depo- 



sition [2^), where a low-energy surface beta can have 
only part of the ionization it generates collected. The 
subset of these episodes where the gamma escapes inter- 
action with another detector could then give rise to a 
family of single-interaction events within the low-energy 
NR band. However, multiples within the ±2a NR detec- 
tor bands and Er<6 keVnr exhibit a fiat distribution in 
ionization energy, as opposed to the singles distribution, 
which is markedly peaked around the NR centroid (^0.6 
keVee) and well-separated from the ER band. In addi- 
tion to this, when the rates of singles and multiples for 
each detector are examined in this same region, they ex- 
hibit a near-perfect lack of linear dependence (Pearson's 
correlation coefficient r = 0.08), also indicating that the 
possibility of a common physical origin is remote. 

The more exotic alternative of a WIMP origin can 
be assessed by comparison to recent anomalies in other 
searches. Fig. 4 overlaps the best-fit to the NR band 
component from the analysis of all summed detectors 
onto CoGeNT data The conversion to ionization 
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FIG. 3: Best-fit components of tfie null and alternative mod- 
els, overlapped on summed data, projected on Er and Ei. 
Red: ER and SE combined through the Crystal Ball PDF 
[isj . Green: ZC. Blue: NR. Black: sum of components. The 
null model requires a large deviation of the ZC centroid to 
Ei ~ 0.5 keVee, hard to reconcile with adequate Ei calibra- 
tions [l2| and with the mean Ei of ZC events above ~ 5 keVnr, 
a region where their true centroid can be assessed (Fig. 2). 
The separation between ZC and NR populations is noticeable 
for data in the 5-11.9 keVnr analysis region used in 



energy is done as in [H, by using the more reliable ger- 
manium quenching factor measured by CoGeNT [23|. 
Assuming this exponential distribution is the approxi- 
mate response to a WIMP, we generate a CDMS region 
of interest (ROI) in WIMP couphng vs. WIMP mass (Fig. 
4, inset) that includes present uncertainties. 

Our analysis allows for a straightforward estimate of 
the sensitivity of the search for an annual modulation in 
[H, by integrating best- fit signal (NR) and background 
(ER, SE, ZC) components inside CDMS's ±2a NR, 5- 
11.9 keVnr, "signal box" (Fig. 1, blue enclosure). We find 
that out of ~ 167 events within, only 35% would corre- 
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FIG. 4: Blue: best-fit NR component for CDMS summed 
detector data, translated to ionization scale and overlapped 
on histogrammed CoGeNT data [3| after normalization to 
the vertical scale. Neither is corrected for efficiency next to 
threshold. Dotted blue lines represent the la uncertainty in 
the parameter Ai for NR. A dashed black line represents 
known CoGeNT backgrounds (flat-|-cosmogenic which 
provide an adequate fit to the data down to ~ 1.2 keVee, 
the lower boundary of the CDMS annual modulation search 
region. Inset: 90% and 99% C.L. CDMS ROI in WIMP cou- 
pling vs. mass (see text), including all present uncertainties 
except for those related to CDMS's energy scales [23|. ROI's 
for CoGeNT [4], CRESST [H and DAMA 0] are from 0, 
and include the effect of a residual surface event contamina- 
tion in CoGeNT described in [1]. The DAMA ROI assumes 
a Max;wellian dark matter halo: deviations from it can dis- 
place it to lower WIMP couplings [1, i, HI, [H] . Additional 
uncertainties for DAMA exist l2ql. 



spond to the putative WIMP (NR) signal. This trans- 
lates into 0.035 NR events /keVnr kg day, whereas the 
99% exclusion claimed in [l| is for modulations larger 
than 0.06 events / keVnr kg day. In other words, even 
at a modulation amplitude of 100%, the search in [l[ 
would fail to exclude a WIMP origin for the NR excess 
seemingly present in CDMS data. Important additional 
concerns about the search in [l[ can be listed. For in- 
stance, the addition of non-overlapping time periods [27 1 
from detectors spanning an order of magnitude in back- 
ground rate within the signal box IJ, |l6| , the negligible 
overlap with the CoGeNT spectral region containing a 
clear excess of events (Fig. 4), or unresolved issues related 
to CDMS's energy scales [23|. However, we emphasize 
that the choice of signal box boundaries, one that results 
in a poor signal-to-background ratio, is already sufficient 
to cripple its sensitivity. 

In conclusion, we find that recently released 2007-2008 
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CDMS data strongly support {5.7a C.L.) the presence of 
a family of low-energy events concentrated in the nuclear 
recoil band. An origin in neutron scattering is highly un- 
likely. Their rate and spectral shape provides a match 
to the majority of low-energy events unaccounted for by 
CoGeNT (Fig. 4). Both searches employ the same tar- 
get material, germanium, but perform orthogonal back- 
ground cuts [8|, enhancing the possible meaning of this 
coincidence. If this excess is interpreted as a WIMP sig- 
nal, it falls in close proximity to other anomalies reported 
by recent dark matter experiments (Fig. 4). The favored 
region of WIMP mass is also of present interest in indi- 
rect searches for dark matter [13 . We also determine that 
the recent search for an annual modulation signal by the 
CDMS collaboration is insufficiently sensitive to exclude 
a dark matter origin for this excess, due to an inadequate 
selection of analysis re gion . Unsupported quantitative 
statements made in 0,]ll| about background composi- 
tion in CDMS detectors are not compatible with our 
findings. 

We encourage the CDMS collaboration to continue re- 
fining their understanding of detector backgrounds and 
energy scales. We are indebted to M. Bellis, D. Hooper, 
C. Kelso, D. Moore, A.E. Robinson and N. Weiner for 
useful discussions. N.F. is supported by the NNSA 
Stewardship Science Graduate Fellowship program un- 
der grant number DE-FC52-08NA28752. 
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